Quark mass and isospin dependence of the deconfining critical temperature 
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We describe the deconfining critical temperature dependence on the pion mass and on the isospin 
chemical potential in remarkably good agreement with lattice data. Our framework incorporates 
explicit dependence on quark masses, isospin and baryonic chemical potentials for the case of two 
flavors through ingredients from well-known high- and low-energy theories. In the low-energy sector, 
the system is described by a minimal chiral perturbation theory effective action, corresponding to 
a hot gas of pion quasiparticles and heavy nucleons. For the high-temperature sector we adopt a 
simple extension of the fuzzy bag model. We also briefly discuss the effects of mass asymmetry and 
baryon chemical potential. 



Introduction. The phase diagram of quark matter has 
been the object of intense investigation during the last 
years, and yet several open questions within the ther- 
modynamics of strong interactions still remain unsolved 
[J, Q]. In this quest, Lattice QCD represents the main 
non-perturbative approach within the full theory [3] , al- 
ways complemented by effective models 

In this article we investigate the effects of finite quark 
masses and isospin number on the equation of state of 
hot and dense strongly interacting matter and on the de- 
confining phase transition within a framework inspired 
by chiral perturbation theory (^PT) and lattice results 
for the pressure and the trace anomaly. The setting we 
propose is simple and completely fixed by vacuum QCD 
properties (measured or simulated on the lattice) and 
lattice simulations of finite temperature QCD. More ex- 
plicitly, there is no fitting of mass or isospin chemical 
potential dependence at all. Nevertheless, our findings 
for the behavior of the critical temperature as a func- 
tion of both the pion mass and the isospin chemical po- 
tential are in remarkably good agreement with lattice 
data. It is crucial to note that several detailed studies 
of chiral models failed to describe Tc(to^) d, H, while 
Polyakov-loop models, whose predictions can be fitted to 
the lattice points for Tc(m^) cannot at the present 
form address isospin effects. In this approach, we can 
also investigate in a straightforward manner the effects 
of quark-mass asymmetry and nonzero baryon chemical 
potential, physical cases in which the Sign Problem devel- 
ops, constraining systematic lattice studies. The predic- 
tions within our framework for these regimes which are 
not yet fully probed by lattice QCD are left for a longer 
publication (for preliminary results, see Q). Through- 
out the present paper, the baryon chemical potential is 
fixed to zero. It is not our aim in this short paper to 
provide a complete description of the full richness of the 
QCD phase diagram. Still, we present a description of 
the dependence of the critical deconfining temperature 
on the quark-average mass and on the isospin chemical 
potential simultaneously. 

Small quark masses and a nonzero baryon chemical po- 



tential have always represented a major challenge for lat- 
tice simulations. Presently, although viable lattice sizes 
still prevent extensive and precise studies with realistic 
quark masses, and the Sign Problem considerably re- 
stricts the applicability of Monte Carlo simulations to 
the description of chemically asymmetric media. Lattice 
QCD is starting to provide results with smaller quark 
masses, and pro bing a larger domain at finite chemical 
potential [13, Hi- For the latter, it is believed that 
much can be learned from simulations of realizations of 
QCD that avoid the Sign Problem, such as those with 
vanishing baryon chemical potential and finite isospin 
density, which has a positive fermionic determinant. 

QCD at finite isospin density is certainly, but not only, 
a playground to test numerical approaches to the case of 
finite baryon density on the Lattice It is also part 

of the physical phase diagram for strong interactions, and 
exhibits a very rich phenomenology [14| . It has been un- 
der careful study during the last years, both theoretically 
and experimentally, with a clear identification of certain 
phenomena that depend directly on the isospin asymme- 
try of nuclear matter at intermediate-energy heavy ion 
experiments (see Ref. [l5| and references therein). 

Nevertheless, theoretical and phenomenological stud- 
ies often focus on the chiral limit of QCD, putting aside 
effects from finite quark masses, and on isospin sym- 
metric hot matter, mainly stimulated by the physical 
scenario found in current high-energy heavy ion colli- 
sion experiments [l6| and the quark-gluon plasma 
Some exceptions are, however, chiral model analyses of 
pion-mass dependence of the finite temperature transi- 
tion (e.g., @,H, 01), which fail completely to reproduce 
the lattice behavior, and Polyakov-loop models predic- 
tions as mentioned above. At finite isospin chemical po- 
tential, though there are phase-diagram investigations in 
PN JL models [3l , none of them has addressed the criti- 
cal temperature dependence on the isospin chemical po- 
tential in particular. 

To investigate the effects of nonzero quark masses 
and isospin asymmetry on the deconfining transition, we 
build an effective theory that combines ingredients from 
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xPT in the low-energy sector with the phenomenological 
fuzzy bag model at high energy. The high-energy regime 
is described perturbatively by two-flavor QCD with mas- 
sive quarks and explicit isospin symmetry breaking. Non- 
perturbative (confinement) effects at this scale are incor- 
porated through a fuzzy bag description [l^ with coef- 
ficients extracted from lattice data. For the low-energy 
sector, we adopt an effective action inspired by xPT that 
exhibits the same structure of symmetries contained in 
the high-energy theory [i^l . The quasi-nucleon degrees of 
freedom described in this regime seem to be crucial for 
understanding the pion-mass dependence of the dccon- 
fining temperature. The definition of parameters, as well 
as masses, is such that variations in the deconfined sec- 
tor arc totally consistent with variations in the confined 
one, which guarantees the bookkeeping in the different 
degrees of freedom that are present in the description. 

By matching the two branches of the equations of state, 
corresponding to the high and low temperature regimes, 
we of course obtain a first-order transition. Recently im- 
proved Lattice QCD calculations, with almost realistic 
quark masses, seem to indicate a crossover instead Q. 
On the other hand, from the experimental standpoint a 
weakly first-order transition is not ruled out, and in fact 
corresponds to the scenario adopted in very successful 
hydrodynamic calculations Although this is a cru- 
cial question for the understanding of the phase structure 
of QCD, it is essentially of no consequence to the analy- 
sis we undertake. The value of the critical temperature 
we obtain is ~ 5 — 10% different from current values 
extracted from the Lattice (^(sij, which is always the 
case in effective field theory approaches and of no harm 
to our analysis, either. Our concern is with providing a 
good description of the behavior of the critical tempera- 
ture for increasing values of quark masses and the isospin 
chemical potential. 

Low-energy sector. The physical setting for the low- 
energy regime of strong interactions is that of a system of 
heavy nucleons in the presence of a hot gas of pions whose 
masses are already dressed by corrections from tempera- 
ture, isospin chemical potential and quark masses. The 
effective Lagrangian reads Ceff = Cn + ^tt + ^ttN, where 

Cn = iVU^-Mw + ^M/7o^^ + ^MB7oW, (1) 



(2) 



C^N = ^N^JJ0^^-l-^lno[r^7^]] N, (3) 



decay rate of the neutron, is the pion decay constant, 
and /i is a function of temperature and isospin chemi- 
cal potential. Nucleons are represented by = {p, n) 
with p, n being the proton and neutron spinors, respec- 
tively, and have a mass matrix Mjv = diag(Mp, M„) = 
diag(Af - SM,M + SM). This corresponds to the 
nucleon chiral Lagrangian [2lj , but considering mass cor- 
rections at zero temperature and chemical potential, and 
the coupling to dressed pions. 

The effective (dressed) masses of the pions tt" ~ tt^ 
and TT^ = '^{'^^ T ^i"^); which depend on tempera- 
ture T, isospin chemical potential /i/ and mass asym- 
metry 6m = {md — mu)/^, are denoted, respectively, by 
Too = Tt^o{T, /i/, Sm) and m± = m±{T, ijj) [32| and their 
explicit expressions were calculated in Refs. |20[. In the 
so-called first phase, a regime in which < m^r, they 
have the form 



Too = Tot, 



1 

1 + -QCTOO 



TO^± = m± ^ h = m-jj 



1 1 (To 

1 + -aai ± -a 

2 2 TO^ 



(4) 

T (5) 



up to first order in a = {m-^ / An /t^)'^ . Here, croo: croi 
CTi, and h are functions of temperature, isospin chemical 
potential and quark masses [20| . 

In the second phase (|/x/| > to^t), a condensation of 
pions occurs, and a superfluid phase sets in [s^. In this 
new phase, in order to reestablish the vacuum structure, 
a chiral rotation is produced due to the isospin symmetry 
breaking. All this produces a pion mixing, and the nucle- 
ons also couple in a different way. The degrees of freedom 
do not correspond anymore to pions, but we can still call 
them quasi-pions since their masses in the two phases 
match at the transition point. The tree level masses do 
not have the shape as in the equations above. Instead, 
= l/z/l, m^- = 0, and = ^/^l + 3(TO^//i/)4 

14|. However, it is possible to treat this phase with a 



where qa is the axial vector current coefficient of the nu- 
cleon, which accounts for renormalization in the weak 



simple approximation near the superfluid phase transi- 
tion. In the regime in which \fii\ > to^t, the natural 
expansion parameter is given by = 1 — m'^/fij [20| . 
after scaling all the parameters by The result of the 
first terms in this expansion (s^ = 0) provides the same 
equations as in the normal phase, just replacing by 

Strictly speaking, this is valid only for values of 
very close to to^t, i.e. < \/S>/7 to^ (for a more de- 
tailed discussion cf. [20[). For simplicity, we also apply 
this results for slightly higher values of the isospin chem- 
ical potential [sJl- These results in are confirmed 
by a NJL analysis [l^l . 

The direct effect of the baryonic chemical potential in 
the pure pion quasiparticle gas is omitted since, without 
considering gluonic corrections, it appears as an anoma- 
lous term in the 0{P'^) chiral Lagrangian [2^ , and will be 
present only in two-loop corrections according to power 
counting. For very large values of /is, one has in princi- 
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pie to incorporate effects from the color superconductiv- 
ity gap in the calculation of meson masses in an effective 
theory near the Fermi surface 
analysis, we treat the case = 



2J, |2g. In the present 



0. 



The nucleon masses, Mp = M — SM and Mn = M + 
SM, are dressed by leading-order contributions in zero- 
temperature baryon xPT. Using the results from Ref. 
[2^, for the isospin symmetric case with explicit chiral 
symmetry breaking, and Ref. 



27| , which includes explicit 



isospin breaking effects, we have (neglecting terms ~ m^, 
^ TTig log(mq), and of higher order in mq): 

M (m) = Mo + 2 71 ?7i + 2^/^ 73/2 m^/^ , (6) 



dM{5m) = 2 7i 



asymm 



5m , 



(7) 



Mq being the nucleon mass in the chiral limit, m = {rriu- 
m(i)/2 the average quark mass, and (ssj 
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(8) 
(9) 
(10) 



Here, all parameters and coefficients are fixed to repro- 
duce properties of the QCD vacuum either measured 
or extracted from recent lattice simulations. Explicitly, 
(qq) = —(225 MeV)fis the (1-flavor) chiral condensate 
in the chiral limit [2^ and, from Table I in Ref. [2^, 
Mo = (0.882 ± 0.003) GeV, ci = (-0.93 ± 0.04) GeY-\ 
gA = 1.267, and U = 92.4 MeV, so that 71 = 4.9630 ± 
0.2135, and 73/2 = -0.273424 MeV^^^ Finally, from 
Table 3 in Ref. [l^l, one can extract {2a — /3)/3. Con- 
verting from lattice units to GeV (1 (lattice units) = 
b = 0.125 fm), in the case of the 0{mq) fit, we arrive 
at 7^^'""^ = 0.16734 ±0.07858. This fixes the dispersion 
relation satisfied by the proton and neutron as 
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Ep/nip) - Vp' + iMT 5MY + -JIB ± , (11) 

where p is the 3-momentum and the antiparticle disper- 
sion relations arc obtained from the ones above by the 
substitution ^ —yn- 

High- energy sector. The fuzzy bag model has been 
proposed by Pisarski [l^ as a phenomenological pa- 
rameterization of the equation of state to account for 
the plateau in the trace anomaly normalized by T^, 
(e — 3p)/r^, observed in lattice results above the crit- 
ical temperature. Besides the usual MIT-type bag con- 
stant, the total pressure for QCD in this model has also 
a non-perturbative contribution ~ [3^: pdccoui{T) ~ 
PpQCD (T) — iJfuzzy T'^ — Bmyt- Thc trace anomaly associ- 
ated with this equation of state, assuming that PpQCD ~ 
is then e - 3p = 2Bfu,,y + 4Bmit. 



Recently, a similar parameterization has been used to 
fit lattice results for thc trace anomaly at high tempera- 
tures 0, T > 1.5 Tc « 300 MeV , yielding 



3p 



2^4 



high T 



= i^05^ + ^ + ^, (12) 



with coefficients 6 and c given in Table VIII of Ref. [9(] . 
Notice that the first term in Eq. p2|) comes from a 0{al) 
perturbative contribution to the pressure and is impor- 
tant for the fit only at very high temperatures. Hence, we 
neglect this term in what follows, obtaining thc following 
values for the bag cocfhcicnts: -Bfuzzy = 0.05 GcV^ and 
Smit = 0.006 GeV. 

In our effective theory, we adopt, phcnomcnologically, 
a simple extension of the fuzzy bag equation of state 
which includes the infiuence of finite chemical potentials 
and masses in the perturbative pressure, neglecting for 
simplicity non-perturbative contributions due to the fi- 
nite quark chemical potentials /i/, so that: Pdcconf — 

PpQCbIT, ^J.f,mf) — Bfuzzy T'^ - Bmit- 

Results for the critical temperature As detailed above, 
we constructed a model using results of well-studied the- 
ories and choosing carefully the relevant ingredients to 
study mass and isospin number effects. Both low- and 
high-energy sectors are completely fixed and now we 
turn to the determination of the prediction of this model 
for the behavior of the deconfining critical temperature. 
From our results for the massive free gas contribution of 
the pQCD pressure in the fuzzy bag model at finite tem- 
perature, isospin and baryon number, and the free gas 
pressure of quasi-pions and nucleons in the low-energy 
regime, the critical temperature and chemical potential 
for the deconfining phase transition are extracted by 
maximizing the total pressure. The validity of our ap- 
proach is, of course, restricted by the scale of xPT: e.g. 
for « rrip « 770 MeV, the expansion parameter in 
xPT is a « 0.45, so that the extension of thc predictions 
to < 1 GeV is justified [stI ]. 

The pion mass dependence, or equivalently the quark 
average mass dependence, of the critical temperature is 
displayed in Fig. [l]for vanishing chemical potentials with 
the temperature and the pion mass normalized by the 
square root of the string tension. Our curves stop at a 
point from which our xPT approach clearly breaks down. 
The results of our framework are compared to lattice data 
from Refs. ^ {Np = 2,3) and ^ {Np = 2 + 1) and 
two other phenomenological treatments: the 0(4) linear 
sigma model 5| and a renormalization group improved 
computation 6| (cf. also Ref. that discusses the 
quark-meson model using the proper-time renormaliza- 
tion group approach). The approximate mass indepen- 
dence observed in the lattice data is very well reproduced 
within our framework, while the other descriptions tend 
to generate a qualitatively different behavior. This fea- 
ture is yet another indication that the functional depen- 
dence of Tc(m7r) requires confinement ingredients to be 
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FIG. 1: Critical temperature as a function of normal- 
ized by the string tension = 425 MeV. Our results are 
compared to lattice data [gl. |29|] and other approaches [1, 0|. 



reproduced, being incompatible with a phase transition 
dictated by pure chiral dynamics. This argument goes in 
the same direction of Ref . Q , the main difference here 
being the fact that we construct the mass dependence 
from the niquark = hmit, with the heavy quasi- nucleons 
as the key new element at low-energies. Moreover, our 
results are not strongly sensitive to the choice between 
the fuzzy bag model and the usual MIT bag model. The 
critical values for the MIT bag model are systematically 
lowered, but the qualitative behavior is not altered, as 
illustrated in Fig. [TJ This indicates that a consistent 
treatment of the quark mass dependence connecting both 
perturbative regimes of energy is probably the essential 
ingredient to describe this observable. 
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FIG. 2: Critical temperature as a function of ^/ for two dif- 
ferent values of the pion mass. Our results are compared to 
lattice data and other approaches [30l |. 

In Fig. O the critical temperature is plotted as a 
function of the isospin chemical potential. The critical 



temperature is normalized by its value in the absence of 
chemical asymmetry, whereas the isospin chemical poten- 
tial is normalized by the critical temperature itself. The 
curves obtained within our framework are represented by 
solid lines. Once again our results are in good agreement 
with lattice computations even though the curve 

that is closer to the lattice points corresponds to a small 
vacuum pion mass, which is not the situation simulated 
on the lattice [3]. Our curves for different values of m.^- 
start to disagree more appreciably for [ii > 2Tc. Re- 
call that our treatment is valid up to isospin chemical 
potentials that are larger but not much larger than the 
pion mass, and contains the effects from pion condensa- 
tion for /i/ > TOtt- Fig- m also exhibits, for comparison, 
results using the hadron resonance gas model [30] that 
appear to depart from the lattice data at a much lower 
value of the isospin chemical potential. The two curves 
are produced by two different methods to determine the 
critical temperature: the dotted curve is obtained from 
the observation that the deconfined phase emerges at a 
constant energy density, whereas the dashed one uses the 
fact that the quark-antiquark condensate for the light 
quarks almost disappears at the quark-hadron transition 
|30|. Similarly to what we observe for the quark-mass 
dependence of Tc, it is clear from Fig. [2] that plain chi- 
ral considerations render the largest discrepancies for the 
behavior of Tdp-i) as compared to lattice data. 

A up-down quark mass imbalance, characterized by the 
relative difference in quark mass, which is in principle 
not much smaller than one, tends to increase the critical 
temperature, though by a quantitatively small amount, 
as expected. The value of the critical baryonic chemical 
potential, beyond which matter is deconfined, also seems 
to increase with the vacuum pion mass. Detailed results 
in these and other thermodynamic observables will be 
presented in a longer pubhcation (see also 3). 

Final remarks. We constructed a frugal effective 
framework selecting features from established theories 
that are relevant in the determination of the quark-mass 
and chemical-asymmetry dependence of the deconfining 
transition. From a simple free gas calculation of the 
equation of state, we found surprisingly good agreement 
with different lattice data for the behavior of the critical 
temperature for the deconfining transition with masses 
and the isospin chemical potential, indicating that the 
model captures the essential features brought about by 
the inclusion of those effects. It should be emphasized 
that the framework provides simultaneous predictions for 
these two functional dependences of the critical tempera- 
ture, TcirriT^) and Tc{fJ,i), with only one set of coefficients 
completely determined by observed QCD vacuum prop- 
erties and lattice simulations at finite temperature and 
zero chemical potentials and fixed quark masses. As far 
as we are aware of, this is the first framework to well- 
reproduce lattice data. Our description of Tc(m.n.) rein- 
forces, through a completely different approach, the dis- 
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cussion in Rcf. Q in which the quark-mass dependence of 
the finite temperature transition on the lattice is shown 
to be compatible with results from slightly perturbed 
Polyakov-loop models, contrasting with the failure of dif- 
ferent chiral models. Confinement properties seem to in- 
fluence strongly the mass-dependence of the QCD transi- 
tions observed on the lattice. Moreover, the approximate 
framework built in the present paper, constituted of se- 
lected ingredients from well-kown high- and low-energy 
theories, is able to access physical settings including 
quark-mass asymmetry and finite baryon chemical poten- 
tial, regimes in which lattice simulations encounter severe 
difficulties due to non-positive definite fermion determi- 
nants. Therefore, this model provides a pragmatic tool 
to investigate the role played by nonzero quark masses 
and chemical potentials, also going beyond the free gas 
approximation, complementing in a healthy direction re- 
sults from other model approaches and lattice simula- 
tions. Since the (mass-symmetric) isospin chemical po- 
tential regime does not suffer from the Sign problem, the 
road is open for detailed Monte Carlo studies and further 
comparisons. On the experimental side, intermediate- 
energy experiments in nuclear physics are providing data 
and observables that are sensitive to chemical asymme- 
try, whereas high-energy heavy ion collisions to start soon 
at RHIC-BES and FAIR-GSI will probe a region of the 
phase diagram of QCD where effects from /is become im- 
portant. Here we have focused on the influence of finite 
quark masses and isospin chemical potential. Predictions 
of this framework for the flnite regime will be reported 
in the near future. This work was partially supported by 
ANPCyT, CAPES, CNPq, FAPERJ, and FUJB/UFRJ. 
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